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Abstract 

This paper considers the classical error correcting problem which is frequently dis- 
cussed in coding theory. We wish to recover an input vector / 6 R™ from corrupted 
measurements y = Af + e. Here, A is an m by n (coding) matrix and e is an arbitrary 
and unknown vector of errors. Is it possible to recover / exactly from the data yl 

We prove that under suitable conditions on the coding matrix A, the input / is the 
unique solution to the £i-minimization problem (||x||f 1 := \xt\) 

min \\y-Ag\\ tl 

provided that the support of the vector of errors is not too large, ||e||f := \{i : e% ^ 
0}| < p ■ m for some p > 0. In short, / can be recovered exactly by solving a simple 
convex optimization problem (which one can recast as a linear program). In addition, 
numerical experiments suggest that this recovery procedure works unreasonably well; 
/ is recovered exactly even in situations where a significant fraction of the output is 
corrupted. 

This work is related to the problem of finding sparse solutions to vastly underde- 
termined systems of linear equations. There are also significant connections with the 
problem of recovering signals from highly incomplete measurements. In fact, the results 
introduced in this paper improve on our earlier work [5]. Finally, underlying the suc- 
cess of l\ is a crucial property we call the uniform uncertainty principle that we shall 
describe in detail. 

Keywords. Linear codes, decoding of (random) linear codes, sparse solutions to under- 
determined systems, l\ minimization, basis pursuit, duality in optimization, linear pro- 
gramming, restricted orthonormality, principal angles, Gaussian random matrices, singular 
values of random matrices. 
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1 Introduction 



1.1 Decoding of linear codes 

This paper considers the model problem of recovering an input vector / G R n from cor- 
rupted measurements y = Af + e. Here, A is an m by n matrix (we will assume throughout 
the paper that m > n), and e is an arbitrary and unknown vector of errors. The problem 
we consider is whether it is possible to recover / exactly from the data y. And if so, how? 

In its abstract form, our problem is of course equivalent to the classical error correcting 
problem which arises in coding theory as we may think of A as a linear code; a linear code 
is a given collection of codewords which are vectors a\, . . . ,a n 6 R m — the columns of the 
matrix A. Given a vector / G R ra (the "plaintext") we can then generate a vector Af in 
R m (the "ciphertext" ) ; if A has full rank, then one can clearly recover the plaintext / from 
the ciphertext Af. But now we suppose that the ciphertext Af is corrupted by an arbitrary 
vector e G R m giving rise to the corrupted ciphertext Af + e. The question is then: given 
the coding matrix A and Af + e, can one recover / exactly? 

As is well-known, if the fraction of the corrupted entries is too large, then of course we have 
no hope of reconstructing / from Af + e; for instance, assume that m = 2n and consider 
two distinct plaintexts /, /' and form a vector g G R m by setting half of its m coefficients 
equal to those of Af and half of those equal to those of Af. Then g = Af + e = Af + e! 
where both e and ef are supported on sets of size at most n = m/2. This simple example 
shows that accurate decoding is impossible when the size of the support of the error vector 
is greater or equal to a half of that of the output Af. Therefore, a common assumption in 
the literature is to assume that only a small fraction of the entries are actually damaged 

\\e\\ £o := \{i : e t ^ 0}\ < p ■ m. (1.1) 

For which values of p can we hope to reconstruct e with practical algorithms? That is, with 
algorithms whose complexity is at most polynomial in the length m of the code A? 

To reconstruct /, note that it is obviously sufficient to reconstruct the vector e since knowl- 
edge of Af + e together with e gives Af, and consequently / since A has full rank. Our 
approach is then as follows. We construct a matrix which annihilates the m x n matrix A 
on the left, i.e. such that FA = 0. This can be done in an obvious fashion by taking a 
matrix F whose kernel is the range of A in R m , which is an n-dimensional subspace (e.g. F 
could be the orthogonal projection onto the cokernel of A). We then apply F to the output 
y = Af + e and obtain 

y = F(Af + e) = Fe (1.2) 

since FA = 0. Therefore, the decoding problem is reduced to that of reconstructing a 
sparse vector e from the observations Fe (by sparse, we mean that only a fraction of the 
entries of e are nonzero). 
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1.2 Sparse solutions to underdetermined systems 

Finding sparse solutions to underdetermined systems of linear equations is in general NP- 
hard [11, 27]. For example, the sparsest solution is given by 

(P ) min \\d\\e subject to Fd = y (= Pe), (1.3) 

ci£R, m 

and to the best of our knowledge, solving this problem essentially require exhaustive searches 
over all subsets of columns of F, a procedure which clearly is combinatorial in nature and 
has exponential complexity. 

This computational intractability has recently led researchers to develop alternatives to 
(Po), and a frequently discussed approach considers a similar program in the £i-norm which 
goes by the name of Basis Pursuit [8]: 

(P) min \\d\\ tl , Fd = y, (1.4) 

where we recall that \\d\\e 1 = YliLi 1^*1- Unlike the ^o-norm which enumerates the nonzero 
coordinates, the ^i-norm is convex. It is also well-known [2] that (Pi) can be recast as a 
linear program (LP). 

Motivated by the problem of finding sparse decompositions of special signals in the field of 
mathematical signal processing and following upon the ground breaking work of Donoho 
and Huo [13], a series of beautiful articles [14, 15, 21, 30] showed exact equivalence between 
the two programs (Po) and (Pi). In a nutshell, this work shows that for m/2 by m matrices 
F obtained by concatenation of two orthonormal bases, the solution to both (Po) and (Pi) 
are unique and identical provided that in the most favorable case, the vector e has at most 
.914-^/ m/2 nonzero entries. This is of little practical use here since we are interested in 
procedures that might recover a signal when a constant fraction of the output is unreliable. 

Using very different ideas and together with Romberg [4], the authors proved that the 
equivalence holds with overwhelming probability for various types of random matrices pro- 
vided that provided that the number of nonzero entries in the vector e be of the order of 
m/logm [5, 6]. In the special case where F is an m/2 by m random matrix with indepen- 
dent standard normal entries, [11] proved that the number of nonzero entries may be as 
large as p-m, where p > is some very small and unspecified positive constant independent 
of m. 

1.3 Innovations 

This paper introduces the concept of a restrictedly almost orthonormal system — a collec- 
tion of vectors which behaves like an almost orthonormal system but only for sparse linear 
combinations. Thinking about these vectors as the columns of the matrix F, we show that 
this condition allows for the exact reconstruction of sparse linear combination of these vec- 
tors, i.e. e. Our results are significantly different than those mentioned above as they are 
deterministic and do not involve any kind of randomization, although they can of course 
be specialized to random matrices. For instance, we shall see that a Gaussian matrix 
with independent entries sampled from the standard normal distribution is restrictedly al- 
most orthonormal with overwhelming probability, and that minimizing the ^i-norm recovers 
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sparse decompositions with a number of nonzero entries of size po ■ m; we shall actually give 
numerical values for p$. 

We presented the connection with sparse solutions to underdetermined systems of linear 
equations merely for pedagogical reasons. There is a more direct approach. To recover / 
from corrupted data y = Af + e, we consider solving the following ^-minimization problem 

(P{) mm \\y-Ag\\ ei . (1.5) 

geti" 

Now / is the unique solution of (P{) if and only if e is the unique solution of (-Pi). In other 
words, (Pi) and (P[) are equivalent programs. To see why these is true, observe on the one 
hand that since y = Af + e, we may decompose g as g = / + h so that 

(Pi) 44> min lie - Ahh, . 

On the other hand, the constraint Fx = Fe means that x = e — Ah for some h G R n and, 
therefore, 

(Pi) 44> min llxlk , x = e — Ah 



44> min lie — Ah\ 



which proves the claim. 



The program (P[) may also be re-expressed as an LP — hence the title of this paper. Indeed, 
the ^-minimization problem is equivalent to 

minl T t, -t<y-Ag<t, (1.6) 

where the optimization variables are t G R m and g G R n (as is standard, the generalized 
vector inequality x < y means that x,- L < y, L for all i). As a result, (P[) is an LP with 
inequality constraints and can be solved efficiently using standard optimization algorithms, 
see [3]. 



1.4 Restricted isometries 

In the remainder of this paper, it will be convenient to use some linear algebra notations. 
We denote by (vj)j^j € BP the columns of the matrix F and by H the Hilbert space 
spanned by these vectors. Further, for any T C J, we let Ft be the submatrix with column 
indices j G T so that 

F T c = ^2c j v j G H. 

To introduce the notion of almost orthonormal system, we first observe that if the columns 
of F are sufficiently "degenerate," the recovery problem cannot be solved. In particular, if 
there exists a non-trivial sparse linear combination Y^jeT c j v j = of the Vj which sums to 
zero, and T = T\ U T2 is any partition of T into two disjoint sets, then the vector y 

jen j<=T 2 
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has two distinct sparse representations. On the other hand, linear dependencies SjeJ CjVj = 
which involve a large number of nonzero coefficients Cj , as opposed to a sparse set of co- 
efficients, do not present an obvious obstruction to sparse recovery. At the other extreme, 
if the (vj)j e j are an orthonormal system, then the recovery problem is easily solved by 
setting Cj = (f,Vj) H - 

The main result of this paper is that if we impose a "restricted orthonormality hypothesis," 
which is far weaker than assuming orthonormality, then (Pi) solves the recovery problem, 
even if the (vj)j^j are highly linearly dependent (for instance, it is possible for m := \ J\ to 
be much larger than the dimension of the span of the Vj's). To make this quantitative we 
introduce the following definition. 

Definition 1.1 (Restricted isometry constants) Let F be the matrix with the finite 
collection of vectors (vj)j^j G R p as columns. For every integer 1 < 5 < \ J\, we define the 
S -restricted isometry constants 5s to be the smallest quantity such that Ft obeys 

(l-<5s)||c|| 2 <||F T c|| 2 <(l + ^)||c|| 2 (1.7) 

for all subsets T C J of cardinality at most S, and all real coefficients (cj)j e x- Similarly, 
we define the 5, S' -restricted orthogonality constants 6s,s' f or S + S' < \J\ to be the smallest 
quantity such that 

\(F T c,F T ,c')\<6s,s>-\\c\\\\c'\\ (1.8) 
holds for all disjoint sets T, T' C J of cardinality \T\ < S and \T'\ < S' . 

The numbers 5s and 9s measure how close the vectors Vj are to behaving like an orthonormal 
system, but only when restricting attention to sparse linear combinations involving no more 
than S vectors. These numbers are clearly non-decreasing in 5, S'. For 5 = 1, the value 
5i only conveys magnitude information about the vectors Vj\ indeed 5\ is the best constant 
such that 

1 - Si < \\vj\\ 2 H < 1 + <5i for all j G J. (1.9) 

In particular, 5± = if and only if all of the Vj have unit length. Section 2.3 establishes 
that the higher 5s control the orthogonality numbers 6s,S'- 

Lemma 1.2 We have 9s,s' < $S+S' < 9s,S' + max(Js, Sg') for all S, S' . 

To see the relevance of the restricted isometry numbers 5s to the sparse recovery problem, 
consider the following simple observation: 

Lemma 1.3 Suppose that S > 1 is such that c^s < 1, and let T C J be such that \T\ < 
S. Let f := Ft c for some arbitrary \T\- dimensional vector c. Then the set T and the 
coefficients (cj)j S T can be reconstructed uniquely from knowledge of the vector f and the 
Vj 's. 

Proof We prove that there is a unique c with ||c||^ < S and obeying / = ^ ■ CjVj. Suppose 
for contradiction that / had two distinct sparse representations f = Ft c = Ft' d where 
\T\,\T'\ < S. Then 

F T ut' d = 0, dj := Cjlj^T - c j'^jeT'- 
Taking norms of both sides and applying (1.7) and the hypothesis 62s < 1 we conclude that 
||g?|| 2 = 0, contradicting the hypothesis that the two representations were distinct. ■ 
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1.5 Main results 



Note that the previous lemma is an abstract existence argument which shows what might 
theoretically be possible, but does not supply any efficient algorithm to recover T and cj 
from / and (vj)j e j other than by brute force search — as discussed earlier. In contrast, 
our main theorem result that, by imposing slightly stronger conditions on 82s, the £±- 
minimization program (Pi) recovers / exactly. 

Theorem 1.4 Suppose that S > 1 is such that 

5 s + s + e s ,2S<l, (1.10) 

and let c be a real vector supported on a set T C J obeying \T\ < S. Put f := Fc. Then c 
is the unique minimizer to 

(Pi) min||d||* Fd = f. 

Note from Lemma 1.2 that (1.10) implies 82s < lj an d is in turn implied by 5s + 82s + £35 < 
1/4. Thus the condition (1.10) is roughly "three times as strict" as the condition required 
for Lemma 1.3. 

Theorem 1.4 is inspired by our previous work [5], see also [4,6], but unlike those earlier 
results, our results here are deterministic, and thus do not have a non-zero probability of 
failure, provided of course one can ensure that the system (vj)j^j verifies the condition 
(1.10). By virtue of the previous discussion, we have the companion result: 

Theorem 1.5 Suppose F is such that FA = and let S > 1 be a number obeying the 
hypothesis of Theorem 1.4- Set y = Af + e, where e is a real vector supported on a set of 
size at most S. Then f is the unique minimizer to 

(P[) mm \\y - Ag\\ ei . 

geB. n 



1.6 Gaussian random matrices 

An important question is then to find matrices with good restricted isometry constants, 
i.e. such that (1.10) holds for large values of S. Indeed, such matrices will tolerate a larger 
fraction of output in error while still allowing exact recovery of the original input by linear 
programming. How to construct such matrices might be delicate. In section 3, however, 
we will argue that generic matrices, namely samples from the Gaussian unitary ensemble 
obey (1.10) for relatively large values of S. 

Theorem 1.6 Assume p < m and let F be a p by m matrix whose entries are i.i.d. 
Gaussian with mean zero and variance 1/p. Then the condition of Theorem 1.4 holds with 
overwhelming probability provided that r = S/m is small enough so that 

r < r*(p, m) 

where r*(p,m) is given in Section 3.5. (By "with overwhelming probability," we mean with 
probability decaying exponentially in m.) In the limit of large samples, r* only depends upon 
the ratio, and numerical evaluations show that the condition holds for r < 3.6 • 10~ 4 in the 
case where p/m = 3/4, r < 3.2 • 1CP 4 when p/m = 2/3, and r < 2.3 • 1CP 4 when p/m = 1/2. 
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In other words, Gaussian matrices are a class of matrices for which one can solve an un- 
derdetermined systems of linear equations by minimizing t\ provided, of course, the input 
vector has fewer than p ■ m nonzero entries with p > 0. We mentioned earlier that this 
result is similar to [11] . What is new here is that by using a very different machinery, one 
can obtain explicit numerical values which were not available before. 

In the context of error correcting, the consequence is that a fraction of the output may 
be corrupted by arbitrary errors and yet, solving a convex problem would still recover / 
exactly — a rather unexpected feat. 

Corollary 1.7 Suppose A is an n by m Gaussian matrix and set p = m — n. Under the 
hypotheses of Theorem 1.6, the solution to (P{) is unique and equal to f. 

This is an immediate consequence of Theorem 1.6. The only thing we need to argue is 
why we may think of the annihilator F (such that FA = 0) as a matrix with independent 
Gaussian entries. Observe that the range of A is a random space of dimension n embedded 
in R m so that the data y = Fe is the projection of e on a random space of dimension p. 
The range of a p by m matrix with independent Gaussian entries precisely is a random 
subspace of dimension p, which justifies the claim. 

We would like to point out that the numerical bounds we derived in this paper are overly 
pessimistic. We are confident that finer arguments and perhaps new ideas will allow to 
derive versions of Theorem 1.6 with better bounds. The discussion section will enumerate 
several possibilities for improvement. 

1.7 Organization of the paper 

The paper is organized as follows. Section 2 proves our main claim, namely, Theorem 1.4 
(and hence Theorem 1.5) while Section 3 introduces elements from random matrix theory 
to establish Theorem 1.6. In Section 4, we present numerical experiments which suggest 
that in practice, (P[) works unreasonably well and recovers the / exactly from y = Af + e 
provided that the fraction of the corrupted entries be less than about 17% in the case 
where m = 2n and less than about 34% in the case where m = An. Section 5 explores the 
consequences of our results for the recovery of signals from highly incomplete data and ties 
our findings with some of our earlier work. Finally, we conclude with a short discussion 
section whose main purpose is to outline areas for improvement. 

2 Proof of Main Results 

Our main result, namely, Theorem 1.4 is proved by duality. As we will see in section 2.2, c 
is the unique minimizer if the matrix Ft has full rank and if one can find a vector w with 
the two properties 

(i) (w,Vj) H = sgn(cj) for all j G T, 

(ii) and \(w, Vj)n\ < 1 for all j ^ T, 
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where sgn(cj) is the sign of Cj (sgn(cj) = for Cj = 0). The two conditions above say that 
a specific dual program is feasible and is called the exact reconstruction property in [5], see 
also [4]. For \T\ < S with S obeying the hypothesis of Theorem 1.4, Ft has full rank since 
5s < 1 and thus, the proof simply consists in constructing a dual vector w; this is the object 
of the next section. 

2.1 Exact reconstruction property 

We now examine the sparse reconstruction property and begin with coefficients {w,Vj)n 
for j T being only small in an £2 sense. 

Lemma 2.1 (Dual sparse reconstruction property, £2 version) LetS,S f > 1 be such 
that 5s < 1, and c be a real vector supported onT C J such that \T\ < S. Then there exists 
a vector w G H such that {w,Vj)n = Cj for all j £ T. Furthermore, there is an "exceptional 
set" E C J which is disjoint from T, of size at most 

\E\ < S', (2.1) 

and with the properties 

I <«;,«,-) I < 9S,S ' 7= ■ \\4 f° r alljgTUE 

and 

j&E S 
In addition, \\w\\jj < K ■ \\c\\ for some constant K > only depending upon 5s- 

Proof Recall that Ft : —>■ H is the linear transformation Ft ct ■= J2j£T c j v j where 

ct '■= (cj)j£T (we use the subscript T in ct to emphasize that the input is a |T|-dimensional 
vector), and let F T be the adjoint transformation 

F T w := ({w,Vj) H )j6T- 

Property (1.7) gives 

1-5 S < X miQ (F T F T ) < X max (F T F T ) <1 + 5 S , 

where A m i n and A max are the minimum and maximum eigenvalues of the positive-definite 
operator F t Ft- In particular, since 5\t\ < 1, we see that F t Ft is invertible with 

IK^r)" 1 !! < ( 2 - 2 ) 
Also note that || ^(i^i^r)- 1 1| < ^1 + 5 S /(1 - 5 S ) and set w £ H to be the vector 

w := Ft^Ft) -1 ^; 



8 



it is then clear that F^w = ct, i.e. (w,Vj)H = Cj for all j G T. In addition, \\w\\ < K ■ \\ct\\ 
with K = \/l + 5s/(l — Ss)- Finally, if T' is any set in J disjoint from T with \T'\ < S' 
and dx' = (dj)jeT' is any sequence of real numbers, then (1.8) and (2.2) give 



\(F^,w,d T ')e 2 (T')\ = \{w,F T ,d T >)t 2 (T') 



= (J2mF T )- 1 c T ) j v j ,Y,d j v j ) H 

j£T j£T> 

< &s,S' ■ IK-^t-^t)" 1 ^!! • \\dr'\\ 

< — H ct • \\dT'\\; 
l-ds 



since dr> was arbitrary, we thus see from duality that 

,S' 

5s 



F T ,w\\ h( T>) < 1 _ 5 \\ct\ 



In other words, 

' \Un „,\I 2 W 2 < . 

I -5 S 



V Ku;,^-)! 2 ) 1 / 2 < — ^-llorll whenever T' C J\T and \T'\ < S'. (2.3) 



jer* 

If in particular if we set 

E := {j G J\T : |<^)| > - ^ • ||c T ||}, 

(1 - ds)Vb' 

then \E\ must obey \E\ < S', since otherwise we could contradict (2.3) by taking a subset 
V of E of cardinality S'. The claims now follow. ■ 

We now derive a solution with better control on the sup norm of \(w,Vj)\ outside of T, by 
iterating away the exceptional set E (while keeping the values on T fixed). 

Lemma 2.2 (Dual sparse reconstruction property, version) Let S > 1 be such 
that Ss + 9s,2S < 1> Q- n d c be a real vector supported onT C J obeying \T\ < S. Then there 
exists a vector w G H such that (w, Vj)n = Cj for all j G T. Furthermore, w obeys 

l^'^l ^ n x —a V7^ • II C H f0r allj ? T - ^ 

(l-ds ~ 9s,2S)VS 

Proof We may normalize X^gtI ?! 2 = v^- Write To := T. Using Lemma 2.1, we can 
find a vector w\ G H and a set T\ C J such that 

T n Ti = 
< S 

(wi,Vj) H = cj for all j G T 

I (wi,Vj) H \ < .f 3 ' 3 ' ; for all j?T U T x 
(1 - o s ) 

j'eTi 5 
ll^i ||# < 
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Applying Lemma 2.1 iteratively gives a sequence of vectors w n+ \ G H and sets T n+ i C J 
for all n > 1 with the properties 

r n n(r uT n+1 ) = 

|^Tl+l| < 5 1 

{w n+ i,Vj) H = (w n ,Vj) H for all j G T n 
{w n+1 ,Vj) H = for all j G T 

| K +1 , I < j^j- (j^A for all j T U T n U T n+1 

IK + ik<( T ^) k 

By hypothesis, we have < 1. Thus if we set 



oo 

W 

n=l 



then the series is absolutely convergent and, therefore, w is a well-defined vector in H. We 
now study the coefficients 

oo 

(w,v j ) H = Y i (-l) n - 1 (w n ,v j ) H (2.5) 

n=l 

for j G J. 

Consider first j G To, it follows from the construction that {wi,Vj)u = Cj and (w n , Vj)n = 
for all n > 2, and hence 

(u>, = Cj for all j G To. 

Second, fix j with j G" To and let ij := {n > 1 : j G T n }. Since T ra and T n+ i are disjoint, 
we see that the integers in the set Ij are spaced at least two apart. Now if n G Ij, then by 
definition j G T n and, therefore, 

(w n+1 ,Vj) H = (w n ,Vj) H . 

In other words, the n and n + 1 terms in (2.5) cancel each other out. Thus we have 

(w,Vj) H = (-l) n_1 (w„,^)//. 

n>l;n,n—l£Ij 

On the other hand, if n, n — 1 and n^O, then j T n n T n _i and 
which by the triangle inequality and the geometric series formula gives 

n>l;n,n— 
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In conclusion, 

Os,s 



| (w, Vj ) H - heij (wo,Vj) H \ < - , 

l-ds - Vs,2S 

and since |T| < S, the claim follows. ■ 

Lemma 2.2 actually solves the dual recovery problem. Indeed, our result states that one 
can find a vector w £ H obeying both properties (i) and (ii) stated at the beginning of the 
section. To see why (ii) holds, observe that ||sgn(c)|| = y/\T\ < V~S and, therefore, (2.4) 
gives for all j ^ T 

i / \ i ^ ® s > s . f\T\ ^ ® s > s ^ i 

(l-ds- VS,2S) V A (1 - d S - 0s,2S) 

provided that 5s + 9s s + &s 2S < 1- 



2.2 Proof of Theorem 1.4 

Observe first that standard convex arguments give that there exists at least one minimizer 
d = (dj)j e j to the problem (Pi). We need to prove that d = c. Since c obeys the constraints 
of this problem, d obeys 

\\d\W < \\c\\ tl =5>,-|. (2.6) 

Now take a w obeying properties (i) and (ii) (see the remark following Lemma 2.2). Using 
the fact that the inner product (w, Vj) is equal to the sign of c on T and has absolute value 
strictly less than one on the complement, we then compute 

Nik = I>; + ^-^)l + EKI 

> Y s S n ( c i)( c i + ( d j ~ c i)) + Y d i( w > v o) H 
= Y N + Y^ - c i^ w > V ^ H + Y d i( w ' V ^ H 

= Y n + ( w > Y d 3 v i - Y c ^ 

j£T j£j j<=T 

= Y\c J \ + (w,f-f) 
= Y\°i\- 

Comparing this with (2.6) we see that all the inequalities in the above computation must 
in fact be equality. Since |(«;, Vj)n\ was strictly less than 1 for all j T, this in particular 
forces dj = for all j T. Thus 

Y( d j ~ c j) v j = f ~ f = 0. 

Applying (1.7) (and noting from hypothesis that 5s < 1) we conclude that dj = Cj for all 
j G T. Thus d = c as claimed. This concludes the proof of our theorem. 
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Remark. It is likely that one may push the condition 5s + 9s, s + @S,2S < 1 a little further. 
The key idea is as follows. Each vector w n in the iteration scheme used to prove Lemma 
2.2 was designed to annihilate the influence of w n -\ on the exceptional set T n _i. But 
annihilation is too strong of a goal. It would be just as suitable to design w n to moderate 
the influence of w n -\ enough so that the inner product with elements in T n _i is small 
rather than zero. However, we have not pursued such refinements as the arguments would 
become considerably more complicated than the calculations presented here. 



2.3 Approximate orthogonality 

Lemma 1.2 gives control of the size of the principal angle between subspaces of dimension 
S and S' respectively. This is useful because it allows to guarantee exact reconstruction 
from the knowledge of the 5 numbers only. 

Proof [Proof of Lemma 1.2] We first show that 9s,s> < 5s+S'- By homogeneity it will 
suffice to show that 

l(Ew> Y c 'j' v i')u\ ^ 6 s+s> 
j£T j'eT' 

whenever \T\ < S, \T\ < S', T,T are disjoint, and EjerM 2 = Ej'eT' \ c 'j'? = L Now 
(1.7) gives 

2(1 - 5 S+S <) < II Y c i v J + Y c 'j' v f\\H ^ 2 (! + s s+s>) 

jeT j'eT' 

together with 

2(1 - S s+S ') < ii Y w ~ Y C WH" ^ 2 (! + W), 

jeT j'eT' 

and the claim now follows from the parallelogram identity 

It remains to show that Ss+s' < 9s + 5s- Again by homogeneity, it suffices to establish that 

KEw, Y W>* - 1| < ($s + 9s) 

jet j'ef 

whenever |T| < S + S' and Yljef l c il 2 = 1- ^° P rove this property, we partition T as 
T = TUT' where |T| < 5 and |T'| < S' and write J^jeT | Cj | 2 = q and J2jeT' \ c j\ 2 = 1 ~ a - 
(1.7) together with (1.8) give 

(1 - <5s)a < CjVj, Cj'Vj')H < (1 + *s)a, 

(1 - <S S ,)(1 - a) < cj-uj, ^ Cj/uj-,)^ < (1 + 5 S ')(l - a), 
jeT' j'eT' 

l<E W E W>*l ^ s,s'^ 1/2 (l - «) 1/2 - 
jeT j'eT 
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Hence 

E " ^ ^ ^ a + ^'(1 " «) + 2^« 1/2 (1 - «) 1/2 

jet j'ef 

< max(Ss,Ss') +6 S 

as claimed. (We note that it is possible to optimize this bound a little further but will not 
do so here.) ■ 



3 Gaussian Random Matrices 



In this section, we argue that with overwhelming probability, Gaussian random matrices 
have "good" isometry constants. Consider apbym matrix F whose entries are i.i.d. 
Gaussian with mean zero and variance 1/p and let T be a subset of the columns. We wish 
to study the extremal eigenvalues of F t Ft- Following upon the work of Marchenko and 
Pastur [25], Geman [20] and Silverstein [28] (see also [1]) proved that 

XminfflFr) - (1 - ^7) 2 a.s. 

in the limit where p and \T\ — > oo with 

\T\/p^ 7 < 1. 

In other words, this says that loosely speaking and in the limit of large p, the restricted 
isometry constant 5 (Ft) for a fixed T behaves like 

1 - 6(F T ) < X min (F T F T ) < A ma x (-Ft) < 1 + S(F T ), S(F T ) « 2^7^+ |T|/p. 



Restricted isometry constants must hold for all sets T of cardinality less or equal to S, and 
we shall make use of concentration inequalities to develop such a uniform bound. Note that 
for X" C T, we obviously have 

Amin(-Fr-Fr) < A m i n (Fy,Fy/) < \raax(F T ' Ft 1 ) < A max (i ? j'i ? T) 

and, therefore, attention may be restricted to matrices of size S 1 . Now, there are large devia- 
tion results about the singular values of Ft [29]. For example, letting cr max (Fr) (resp. <r mm ) 
be the largest singular value of Ft so that <r max (-F;r) = A max (i ? yi ? T) (resp. 0- min (-Fr) = 
AminC^r-Pr))) Ledoux [24] applies the concentration inequality for Gaussian measures, and 
for a each fixed t > 0, obtains the deviation bounds 

P (<w(Ft) > 1 + V^ITp + + 1) < e-P t2 / 2 (3.1) 
P (<x min (F T ) < 1 - ^Xr\Tv + o(l) -t)< e'^' 2 - (3.2) 

here, o(l) is a small term tending to zero as p — > oo and which can be calculated explicitly, 
see [16]. For example, this last reference shows that one can select o(l) in (3.1) as —573 ■ 

7 l/6 (1 + v ^ )2 /3. 
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Lemma 3.1 Put r = S/m and set 

f(r) :=^nT V - (^F + y/2H(rfj , 

where H is the entropy function H(q) := — qlogq — (1 — q) log(l — q) defined for < q < 1. 
For each e > 0, the restricted isometry constant 5s of a p by m Gaussian matrix F obeys 

P (1 + 6 S > [1 + (1 + e)f(r)f) < 2 • e -"^W-«/2. (3.3) 

Proof As discussed above, we may restrict our attention to sets \T\ such that |T| = S. 
Denote by r) p the o(l)-term appearing in either (3.1) or (3.2). Put A max = A max (i ? yi ? T) for 
short, and observe that 

P| sup A max > (1 + ^/S/p'+Vp + t) 2 ) < |{T: |T| = 5}|p(A max > (1 + v^ + ^ + t) 2 ) 

\T:\T\=S J V 7 

< Q e -^. 

iFrom Stirling's approximation logm! = mlogm — m + O(logm) we have the well-known 
formula 



771 

log (<-,) = mH(r) + O(logm). 



which gives 



P sup A max > (1 + yfsfv + Vp + tf < e mH ^ ■ e° (l °^ ■ e^ 2 , 

\T:\T\=S J 

The exact same argument applied to the smallest eigenvalues yields 

P ( t m£_ s A min < (1 - 7^ - r? p - t) 2 ) < e m "M • e °( lo sM) • e~^l 2 . 

Fix rjp + t = (1 + e) • y/m/p ■ y/2H(r). Assume now that m and p are large enough so that 
i] p < e/2 • y 7 "^ • V^ H ( r )- Then 

P| sup A max > (1 + 7^ + e- V / 2^)) 2 ) <e" m ^ (r) " /2 - 

\T:|T|=S / 

where we used the fact that the term O(logm) is less than meH(r)/2 for sufficiently large 
m. The same bound holds for the minimum eigenvalues and the claim follows. ■ 

Ignoring the e's, Lemma 3.1 states that with overwhelming probability 

5 s <-l + [l + /(r)] 2 . (3.4) 

A similar conclusion holds for 82s and 83s and, therefore, we established that 

3 

5s + S 2 s + S 3 s <P P /m(r), Pp/m {r) = ^2-l + [l + f{jr)} 2 . (3.5) 

i=i 

with very high probability. In conclusion, Lemma 1.2 shows that the hypothesis of our main 
theorem holds provided that the ratio r = S/m be small so that p p / m (r) < 1. In other 
words, in the limit of large samples, S/m maybe taken as any value obeying p p / m (S/m) < 1 
which we used to give numerical values in Theorem 1.6. Figure 1 graphs the function p p / m ( r ) 
for several values of the ratio p/m. 
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numerical evaluation of the upper bound rho(r) 

3.5 r 

p/m = 3/4 

- - p/m = 2/3 
- p/m = 1/2 

3 - 




1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 

X10- 3 



Figure 1: Behavior of the upper bound p p / m (r) for three values of the ratio p/m, namely, 
p/m = 3/4,2/3, 1/2. 

4 Numerical Experiments 

This section investigates the practical ability of t\ to recover an object / £ R n from 
corrupted data y = Af + e, y € R m (or equivalently to recover the sparse vector of errors 
e € R m from the underdetermined system of equations Fe = z £ R m_n ). The goal here 
is to evaluate empirically the location of the breakpoint as to get an accurate sense of the 
performance one might expect in practice. In order to do this, we performed a series of 
experiments designed as follows: 

1. select n (the size of the input signal) and m so that with the same notations as before, 
A is an n by m matrix; sample A with independent Gaussian entries; 

2. select S percentage of m; 

3. select a support set T of size \T\ = S uniformly at random, and sample a vector e on 
T with independent and identically distributed Gaussian entries 1 ; 

4. make y = Ax + e (the choice of x does not matter as is clear from the discussion and 
here, x is also selected at random), solve (P{) and obtain x*; 

5. compare x to x*; 

6. repeat 100 times for each S and A; 

7. repeat for various sizes of n and m. 

The results are presented in Figure 2 and Figure 3. Figure 2 examines the situation in 
which the length of the code is twice that of the input vector m = 2n, for m = 512 and 

Just as in [6], the results presented here do not seem to depend on the actual distribution used to sample 
the errors. 
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Empirical frequency of exact reconstruction, n = 256, m = 2n 



Empirical frequency of exact reconstruction, n = 512, m = 2n 



0.15 0.2 0.25 0.3 0.35 
# Fraction of corrupted entries 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
# Fraction of corrupted entries 



(a) 



(b) 



Figure 2: Hi-recovery of an input signal from y — Af + e with A anm by n matrix with independent 
Gaussian entries. In this experiment, we 'oversample' the input signal by a factor 2 so that m = 2n. 
(a) Success rate of (Pi) for m = 512. (b) Success rate of (Pi) for m — 1024. Observe the similar 
pattern and cut-off point. In these experiments, exact recovery occurs as long as about 17% or less 
of the entries are corrupted. 



m = 1024. Our experiments show that one recovers the input vector all the time as long 
as the fraction of the corrupted entries is below 17%. This holds for m = 512 (Figure 2(a)) 
and m = 1024 (Figure 2(b)). In Figure 3, we investigate how these results change as the 
length of the codewords increases compared to the length of the input, and examine the 
situation in which m = 4n, with m = 512. Our experiments show that one recovers the 
input vector all the time as long as the fraction of the corrupted entries is below 34%. 



5 Optimal Signal Recovery 



Our recent work [5] developed a set of ideas showing that it is surprisingly possible to 
reconstruct interesting classes of signals accurately from highly incomplete measurements. 
The results in this paper are inspired and improve upon this earlier work and we now 
elaborate on this connection. Suppose we wish to reconstruct an object a in R m from the 
K linear measurements 

yk = (a,4>k) k = l,...,K or y = Fa, (5.1) 

with (frk, the kth. row of the matrix F. Of special interest is the vastly underdetermined 
case, K « N, where there are many more unknowns than observations. We choose to 
formulate the problem abstractly but for concreteness, we might think of a as the coefficients 
a = fy* f of a digital signal or image / in some nice orthobasis, e.g. a wavelet basis so that 
the information about the signal is of the form y = Fa = F$?*f. 

Suppose now that the object of interest is compressible in the sense that the reordered 
entries of a decay like a power-law; concretely, suppose that the entries of a, rearranged in 
decreasing order of magnitude, |a|n) > |a|(2) > • • • > |«| ( m ), obey 

\a\ (k) <B-k- s (5.2) 
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Empirical frequency of exact reconstruction, n = 1 28, m = 4n 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

# Fraction of corrupted entries 

Figure 3: t\-recovery of an input signal from y = Af + e with A anm by n matrix with independent 
Gaussian entries. In this experiment, we 'oversample' the input signal by a factor 4 so that m = An. 
In these experiments, exact recovery occurs as long as about 34% or less of the entries are corrupted. 

for some s > 1. We will denote by J~ S {B) the class of all signals a G R m obeying (5.2). The 
claim is that it is possible to reconstruct compressible signals from only a small number of 
random measurements. 

Theorem 5.1 Let F be the measurement matrix as in (5.1) and consider the solution a" 
to 

min ||a|U, subject to Fa = y. (5.3) 
6eR m 

Let S < K such that 5s + 29s + &S,2S < 1 o,nd set A = K/S. Then a" obeys 

sup ||a-a tt || < C- {K/X)-^- 1 / 2 ^. (5.4) 

aeF s (B) 

To appreciate the content of the theorem, suppose one would have available an oracle letting 
us know which coefficients a^, 1 < k < m, are large (e.g. in the scenario we considered 
earlier, the oracle would tell us which wavelet coefficients of / are large). Then we would 
acquire information about the K largest coefficients and obtain a truncated version ax of 
a obeying 

\\a-a K \\ yc-K-^-W, 

for generic elements taken from J- S {B). Now (5.4) says that not knowing anything about 
the location of the largest coefficients, one can essentially obtain the same approximation 
error by nonadaptive sampling, provided the number of measurements be increased by a 
factor A. The larger S, the smaller the oversampling factor, and hence the connection with 
the decoding problem. Such considerations make clear that Theorem 5.1 supplies a very 
concrete methodology for recovering a compressible object from limited measurements and 
as such, it may have a significant bearing on many fields of science and technology. We 
refer the reader to [5] and [12] for a discussion of its implications. 

Suppose for example that F is a Gaussian random matrix as in Section 3. We will assume 
the same special normalization so that the variance of each individual entry is equal to 1/K. 
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Calculations identical to those from Section 3 give that with overwhelming probability, F 
obeys the hypothesis of the theorem provided that 

S<K/\ X = p ■ log (m/K), 

for some positive constant p > 0. Now consider the statement of the theorem; there is a 
way to invoke linear programming and obtain a reconstruction based upon 0(K log (m/K)) 
measurements only, which is at least as good as that one would achieve by knowing all the 
information about / and selecting the K largest coefficients. In fact, this is an optimal 
statement as (5.4) correctly identifies the minimum number of measurements needed to 
obtain a given precision. In short, it is impossible to obtain a precision of about if _ ( s_1 / 2 ) 
with fewer than K\og(m/K) measurements, see [5, 12]. 

Theorem 5.1 is stronger than our former result, namely, Theorem 1.4 in [5]. To see why 
this is true, recall the former claim: [5] introduced two conditions, the uniform uncertainty 
principle (UUP) and the exact reconstruction principle (ERP). In a nutshell, a random 
matrix F obeys the UUP with oversampling factor A if F obeys 

^ < 1/2, S = p-K/X, (5.5) 

with probability at least 1 — 0(iV _7 / p ) for some fixed positive constant 7 > 0. Second, a 
measurement matrix F obeys the ERP with oversampling factor A if for each fixed subset 
T of size \T\ < S (5.5) and each 'sign' vector c defined on T, there exists with the same 
overwhelmingly large probability a vector w G H with the following properties: 

(i) (w,Vj) = Cj, for all j G T; 

(ii) and \(w,Vj}\ < \ for all j not in T. 

Note that these are the conditions listed at the beginning of section 2 except for the 1/2 
factor on the complement of T. Fix a G J- S (B). [5] argued that if a random matrix obeyed 
the UUP and the ERP both with oversampling factor A, then 

\\a-J\\ < C- {K/Xy (s - 1/2 \ 

with inequality holding with the same probability as before. Against this background, 
several comments are now in order: 

• First, the new statement is more general as it applies to all matrices, not just random 
matrices. 

• Second, whereas our previous statement argued that for each a G R m , one would 
be able — with high probability — to reconstruct a accurately, it did not say anything 
about the worst case error for a fixed measurement matrix F. This is an instance 
where the order of the quantifiers plays a role. Do we need different F's for different 
objects? Theorem 5.1 answers this question unambiguously; the same F will provide 
an optimal reconstruction for all the objects in the class. 

• Third, Theorem 5.1 says that the ERP condition is redundant, and hence the hypoth- 
esis may be easier to check in practice. In addition, eliminating the ERP isolates the 
real reason for success as it ties everything down to the UUP. In short, the ability 
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to recover an object from limited measurements depends on how close F is to an or- 
thonormal system, but only when restricting attention to sparse linear combinations 
of columns. 

We will not prove this theorem as this is a minor modification of that of Theorem 1.4 in 
the aforementioned reference. The key point is to observe that if F obeys the hypothesis of 
our theorem, then by definition F obeys the UUP with probability one, but F also obeys 
the ERP, again with probability one, as this is the content of Lemma 2.2. Hence both the 
UUP and ERP hold and therefore, the conclusion of Theorem 5.1 follows. (The fact that 
the ERP actually holds for all sign vectors of size less than S is the reason why (5.4) holds 
uniformly over all elements taken from T S {B), see [5].) 

6 Discussion 

6.1 Connections with other works 

In our linear programming model, the plaintext and ciphertext had real-valued components. 
Another intensively studied model occurs when the plaintext and ciphertext take values 
in the finite field F<i := {0,1}. In recent work of Feldman et al. [17], [18], [19], linear 
programming methods (based on relaxing the space of codewords to a convex polytope) 
were developed to establish a polynomial-time decoder which can correct a constant fraction 
of errors, and also achieve the information-theoretic capacity of the code. There is thus 
some intriguing parallels between those works and the results in this paper, however there 
appears to be no direct overlap as our methods are restricted to real-valued texts, and the 
work cited above requires texts in F^. Also, our error analysis is deterministic and is thus 
guaranteed to correct arbitrary errors provided that they are sufficiently sparse. 

The ideas presented in this paper may be adapted to recover input vectors taking values 
from a finite alphabet. We hope to report on work in progress in a follow-up paper. 

6.2 Improvements 

There is little doubt that more elaborate arguments will yield versions of Theorem 1.6 with 
tighter bounds. Immediately following the proof of Lemma 2.2, we already remarked that 
one might slightly improve the condition 5 s + &s,S + 6s,2S < 1 at the expense of considerable 
complications. More to the point, we must admit that we used well-established tools from 
Random Matrix Theory and it is likely that more sophisticated ideas might be deployed 
successfully. We now discuss some of these. 

Our main hypothesis reads 5s + 6s,s + #S,2S < 1 but in order to reduce the problem to the 
study of those 5 numbers (and use known results), our analysis actually relied upon the 
more stringent condition 5s + ^25 + 5ss < 1 instead, since 

Ss + &S,S + 6s,2S < $s + fos + $3S- 

This introduces a gap. Consider a fixed set T of size \T\ = S. Using the notations of that 
Section 3, we argued that 

5(F T ) + 
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and developed a large deviation bound to quantify the departure from the right hand-side. 
Now let T and T" be two disjoint sets of respective sizes 5 and 5' and consider 9 (Ft, Ft 1 )'- 
9 (Ft, Ft') is the cosine of the principal angle between the two random subspaces spanned 
by the columns of Ft and Ft' respectively; formally 

9(Ft,Ft>) = sup(u, u'), u G span(Fy), v! G span(Fy), ||u|| = ||«'|| = 1. 

We remark that this quantity plays an important analysis in statistical analysis because of 
its use to test the significance of correlations between two sets of measurements, compare 
the literature on Canonical Correlation Analysis [26]. Among other things, it is known [31] 
that 

9 (Ft, F t >) ->■ Vt(1 - i) + vV(l - l) a.s. 

as p — ► oo with S/p — ► 7 and 5'/p — > 7'. In other words, whereas we used the limiting 
behaviors 

<5(F 2T ) 2^/2^ + 2 7 , <5(F 3T ) -► 2^/3^ + 3 7 , 
there is a chance one might employ instead 

9(F T ,F TI ) ^2^7(1-7), 9(F T ,F TI ) -► V7(l 1 27) + V 2 7(l " 7) 

for |T| = |T'| = 5 and |T'| = 2|T| = 25 respectively, which is better. Just as in Section 3, 
one might then look for concentration inequalities transforming this limiting behavior into 
corresponding large deviation inequalities. We are aware of very recent work of Johnstone 
and his colleagues [23] which might be here of substantial help. 

Finally, tighter large deviation bounds might exist together with more clever strategies to 
derive uniform bounds (valid for all T of size less than 5) from individual bounds (valid for 
a single T). With this in mind, it is interesting to note that our approach hits a limit as 

liminf 5 S + 9 s ,s + 9 S ,2S > J(m/p • r), (6.1) 

S— »oo, S/m^r 

where J(r) := 2^r + r + (2 + y/2)^/r(l - r) + yV(l - 2r). Since J(r) is greater than 1 if 
and only if r > 2.36, one would certainly need new ideas to improve Theorem 1.6 beyond 
cut-off point in the range of about 2%. The lower limit (6.1) is probably not sharp since it 
does not explicitly take into account the ratio between m and p; at best, it might serve as 
an indication of the limiting behavior when the ration p/m is not too small. 

6.3 Other coding matrices 

This paper introduced general results stating that it is possible to correct for errors by 
-£i-minimization. We then explained how the results specialize in the case where the coding 
matrix A is sampled from the Gaussian ensemble. It is clear, however, that one could 
use other matrices and still obtain similar results; namely, that (P[) recovers / exactly 
provided that the number of corrupted entries does not exceed p ■ m. In fact, our previous 
work suggests that partial Fourier matrices would enjoy similar properties [5,6]. Other 
candidates might be the so-called noiselets of Coifman, Geshwind and Meyer [9]. These 
alternative might be of great practical interest because they would come with fast algorithms 
for applying A or A* to an arbitrary vector g and, hence, speed up the computations to 
find the £i-minimizer. 
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